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ABSTRACT 

In the dark matter (DM) halos embedding galaxies and galaxy systems the 'entropy' K = a 2 / p 2 / 3 (a quantity 
that combines the radial velocity dispersion a with the density p) is found from intensive A^-body simulations to 
follow a powerlaw run K oc r a throughout the halos' bulk, with a around 1 .25. Taking up from phenomenology 
just that a w const applies, we cut through the rich analytic contents of the Jeans equation describing the self- 
gravitating equilibria of the DM; we specifically focus on computing and discussing a set of novel physical 
solutions that we name a-profiles, marked by the entropy slope a itself, and by the maximal gravitational 
pull Kcatipi) required for a viable equilibrium to hold. We then use an advanced semianalytic description 
for the cosmological buildup of halos to constrain the values of a to within the narrow range 1 .25 - 1 .29 
from galaxies to galaxy systems; these correspond to halos' current masses in the range 10 11 - 10 15 M Q . Our 
range of a applies since the transition time that - both in our semianalytic description and in state-of-the-art 
numerical simulations - separates two development stages: an early violent collapse that comprises a few 
major mergers and enforces dynamical mixing, followed by smoother mass addition through slow accretion. 
In our range of a we provide a close fit for the relation K cr it(a), and discuss a related physical interpretation 
in terms of incomplete randomization of the infall kinetic energy through dynamical mixing. We also give an 
accurate analytic representation of the a-profiles with parameters derived from the Jeans equation; this provides 
straightforward precision fits to recent detailed data from gravitational lensing in and around massive galaxy 
clusters, and thus replaces the empirical NFW formula relieving the related problems of high concentration and 
old age. We finally stress how our findings and predictions as to a and K, alt contribute to understand hitherto 
unsolved issues concerning the fundamental structure of DM halos. 

Subject headings: Dark matter - galaxies: clusters: general - galaxies: halos 



1. INTRODUCTION 

Galaxies, galaxy groups and clusters are widely held to 
form under the drive of the gravitational instability that acts 
on initial perturbations modulating the cosmic density of the 
dominant cold dark matter (DM) component. The instability 
at first is kept in check by the cosmic expansion, but when the 
local gravity prevails collapse sets in. The standard sequence 
runs as follows: a slightly overdense region expands more 
slowly than its surroundings, progressively detaches from the 
Hubble expansion, halts and turns around; then it collapses, 
and eventually virializes to form a DM 'halo' in equilibrium 
under self-gravity. The amplitude of more massive pertur- 
bations is smaller, so the formation is hierarchical with the 
massive structures forming typically later (see Peebles 1993). 

Such a formation history has been confirmed and resolved 
to a considerable detail by many intensive A^-body simula- 
tions. Early on (White 1986) these added an important block 
of information, i.e., the growth of a halo actually includes 
merging events with other clumps of sizes comparable ('ma- 
jor' mergers) or substantially smaller ('minor' mergers), down 
to nearly smooth accretion (see Springel et al. 2006). 

More recently, a second round of features has been increas- 
ingly recognized in highly resolved simulations of individual 
halos forming in cosmological volumes (Zhao et al. 2003; 
Diemand et al. 2007), to the effect of identifying in the growth 
two definite stages: an early fast collapse including a few vio- 
lent major mergers, and a later calmer stage comprising many 
minor mergers and smooth accretion. During the early col- 
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lapse a substantial mass is gathered, while the major merg- 
ers reshuffle the gravitational potential wells and cause the 
collisionless DM particles to undergo dynamical relaxation; 
therefrom the halo emerges with a definite structure of the in- 
ner density and gravitational potential. During the later stage, 
moderate amounts of mass are slowly accreted mainly onto 
the outskirts, little affecting the inner structure and potential, 
but quiescently rescaling upwards the overall size (see Hoff- 
man et al. 2007). 

In the resulting simulated halos, the coarse-grained proxy 
Q = p/<r 3 to the phase-space density - defined in terms of the 
DM density p and of the radial velocity dispersion a - is em- 
pirically found to be stable after the collapse, and distributed 
in the form of a powerlaw with uniform slope dlogQ/dlogr 
around -1.9 throughout the halo's bulk where it is amenable 
to noise-free measurements (see Taylor & Navarro 2001; 
Dehnen & McLaughlin 2005; Ascasibar & Gottlober 2008). 

Taking up such a notion, the present paper starts (§ 2) from 
investigating the wide set of macroscopic equilibria for the 
DM under self-gravity, as described by the Jeans equation. It 
proceeds (§ 3) to the semianalytic description of their time 
buildup to constrain their governing parameters, and so pre- 
dict the halo realistic structure. It proposes (§ 4) physical 
interpretations for such parameters. The paper is concluded 
(§ 5) with a summary of our findings and a discussion of the 
perspectives they open. 

Throughout we adopt the standard, flat cosmological model 
with normalized matter density Q,m = 0.27, dark energy den- 
sity £1a = 0.73, and Hubble constant Hq = 72 km s" 1 Mpc" 1 . 

2. JEANS EQUILIBRIA 

The static equilibria of the DM halos obey the Jeans equa- 
tion that, on neglecting for the time being anisotropies and the 
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FIG. 1 . — Top panels: the -f(a) and k(o) relations for the powerlaw solutions of the Jeans equation, see § 2.1; blue dots refer to the solution with a = 5/4, and 
green dots refer to the solution with a = 4/3 (the isothermal sphere). In the right panel the dashed line represents the relation ft = 3a/2, and the arrow indicates 
how k increase toward the values (given in the bottom right panel) for the ce-profiles, see § 2.1-2.2 for details. Bottom panels: the 7-7' and K cr j((«) relations for 
the a-profiles with a in the range 1.2 — 35/27, see § 2.1-2.2; blue dots and curves refer to the solution with a = 5/4, and red dots and curves refer to the solution 
with a = 35/27 '. In the left panel, the squares mark the couple of values 70, 7q at r = rQ. Note that we will find in § 3 the relevant values of a to be constrained 
within the range 5/4 — 35/27, i.e., between the blue and red dots/curves. 



blend into the large-scale environment to be discussed in § 5, 
reads 

l d(pa 2 ) _ GM«r) 
p dr r 2 

To elicit the contents of this macroscopic-like equation, 
an explicit link is required between a 2 (r) and p(r). Such a 
link may be provided by Q{r), or by the equivalent quantity 
K(r) = a 2 1 p 2 / 3 = Q~ 2 ^(r); this is often named DM 'entropy' 
in view of its formal analogy with the adiabat of a gas in ther- 
mal equilibrium, and is increasingly used in the context of 
DM equilibria, from Taylor & Navarro (2001) to Faltenbacher 
et al. (2007). We focus on this function of state K(r) in the 
configuration space, and comment on our choice in § 5. 

In terms of the slope dlog,K/dlogr = a, Eq. (1) is recast 
into the particularly simple form 



5 dlog p dlog/f v 2 (r) 
3 dlog r dlog r <J 2 (r) 



(2) 



Here v 2 (r) = GM(< r) jr is the (squared) circular velocity, that 
provides a natural scaling for the depth of the DM potential 
well. 

Taking up from the simulation phenomenology only the no- 
tion that a be uniform throughout the halo bulk (see § 1), 
we will discuss the range of values to be expected for a as 
a clue to their origin. In fact, in their pioneering work Tay- 
lor & Navarro (2001) selected a value of the exponent close 
to a = 1.25. Other, recent studies have provided similar val- 
ues for a: around 1 .27 according to Ascasibar et al. (2004), 
and around 1.29 according to Rasia et al. (2004). Beyond 
the minor differences, it is the narrow empirical range that 
constituted an unsolved challenge. 



2.1. Background 

We begin with looking first (after Taylor & Navarro 2001) 
at powerlaw solutions of Eq. (2) of the form p oc r" 7 , with the 
entropy run set to K oc r a ; then not only a is a given param- 
eter, but also the slope 7 = -dlogp/dlogr is a constant to be 
determined from Eq. (2) in the form 

2-0-7/3 



7 — a = 

' 5 5 3-7 



(3) 



Here the 'bare' constant k = An G pa r^ / cr^ = (3 - 7 )v 2 1 /er 2 1 
compares the gravitational potential to the random kinetic 
energy associated with the velocity dispersion at a refer- 
ence point ro- At the r.h.s. n appears as modulated by 
r-dependencies, produced by the integrated dynamical mass 
M(< r) to yield v c {r) oc r 2_7 /(3-7) in the powerlaw case, and 
by the local dependencies from a 2 (r) oc K(r) p 2 ^(r) oc r a ~ 2j / 3 . 

For a powerlaw solution to actually hold the r.h.s. must be 
constant, which yields the link 

7 = 6-3a, (4) 

as pointed out by Taylor & Navarro (2001), Hansen (2004), 
and Austin et al. (2005). Given this, the balance of the three 
terms in Eq. (3) imposes on the dimensionless gravitational 
pull the constraint 

K = (5 7 /3-a)(3-7) = 6(a-l)(5-3a), (5) 

whose nonlinearities reflect the DM self-gravity. Thus the al- 
lowed powerlaw solutions must lie on a parabola in the plane 
a— k and on a straight line in the plane a -7, as illustrated in 
the top panels of Fig. 1; note that to ensure n > the over- 
all bounds 1 < a < 5/3 must hold, implying that 1 < 7 < 3 
applies for the powerlaw solutions. 
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FIG. 2. — Radial runs for the a-profiles: density p, mass M, circular velocity v c , gravitational potential <E>, velocity dispersion a, and entropy K; the profiles are 
normalized to 1 at the point ro where 7 = 70 = 6 — 3a holds. Various curves are for a in the range 1.2-35/27; the blue one refers to the solution with a = 5/4, 
and the red one refers to a = 35/27. 



Two notable instances correspond to the values a = 4/3 
and 5/4. The former yields the standard, singular isothermal 
sphere with 7 = 2 and n = 2; the latter holds for the value of a 
originally selected by Taylor & Navarro (2001), and implies 
7 = 9/4 and k = 15/8. Note that they are located on the in- 
tersections of the parabola Eq. (5) with the line k = 3 a/2, see 
Fig. 1. 

The powerlaw solutions of the Jeans equation with slope 
70 = 6 - 3 a will turn out to be tangent to the full solutions at 
r = ro so as to provide a useful reference slope in the mid- 
dle of the halo. However, they do not qualify as relevant 
solutions since at both extremes they suffer of unphysical 
features: since 1 < 7 < 3 holds as shown above, the force 
GM(< r)/r 2 oc r 1 " 7 would diverge at the center, while the 
mass M(< r) oc r 3 " 7 would diverge toward large r. 

To avoid these drawbacks, in all physically relevant solu- 
tions the running slope 7(7) is to bend down from the middle 
value in going toward large and small r. In fact, for r <C ro one 
can infer from Eq. (3) that asymptotically 7 — > 7,, = 3 a/5 is to 
hold if the gravitational force is to be finite, which implies that 
both sides of Eq. (3) vanish (see Austin et al. 2005); on the 



other hand, for r ^> ro it is seen that either 7— >7k = 3(l+a)/2 
holds or a density cutoff occurs, if the mass M{< r) is to sat- 
urate to a constant, implying v^(r) oc l/r. In other words, 
powerlaws with inner j a and outer slope 7/, will constitute 
asymptotic or effective approximations to the full solutions. 
For example, for the value a = 5/4 one recovers the slopes 
7„ = 3/4 and j b = 27 /8, found by Taylor & Navarro (2001) in 
their archetypical full solution. 

Thus all full solutions may be conveniently visualized as 
bent down from the corresponding tangent powerlaw so as to 
meet the above asymptotic forms at both small and large radii. 
These solutions will be analytically studied and constrained 
next. In fact, the solution space of the Jeans equation is an- 
alytically rich, so we refer the interested reader to the works 
by Austin et al. (2005), Dehnen & McLaughlin (2005), and 
Barnes et al. (2006, 2007) for extensive coverage. Here we 
summarize the main results and proceed to focus on the 'criti- 
cal' solutions in the terminology of Taylor & Navarro (2001). 

To highlight the properties of these solutions with run- 
ning 7(r), it is convenient (see Dehnen & McLaughlin 2005; 
Austin et al. 2005; also Barnes et al. 2006, 2007) to replace 
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FIG. 3. — Energy distributions for the ct-profiles: specific potential energy £ po t, random kinetic energy £kin> virial ratio 2£| c j„/|£'p t|, and ratio of circular 
velocity to the velocity dispersion v^/a 2 appearing in Eq. (2). Lines as in Fig. 2. 



the original Jeans equation 7 - 3 a/5 = 3 v 2 c /5 a 2 with its first 
and second derivative, and use the differential mass definition 
dM/dr = 4npr 2 , to obtain 

2 3 / t \ f \ 

7'-3(7-7«)(7-7>)=5«(-) [-) , 

(6) 

7 "7 I 7 ^ I = g(7-7a)(7-7o)(7-7W ■ 

Here 7' and 7" denote the first and second logarithmic deriva- 
tive of 7 with respect to r; in addition, ro is the point where 
the density slope matches 70 = 6 -3 a (see Eq. [4]), while 
7 fl = 3 a/5 and 7fe = 3(l + a)/2 are shorthands for the inner 
and the outer effective slopes anticipated above. 

In the second equation the normalization constant n = 
AnGpor^j 'erg does not appear explicitly; however, from the 
first Eq. (6) evaluated at ro it is easily recognized that 

K = 6(a-l)(5-3a) + 5 7 /3 (7) 

is related to 7 , and in a full solution specifically measures the 
upward deviation of k from the corresponding powerlaw, that 
is, above the parabola expressed by Eq. (5), see also Fig. 1 . 

In fact, as the value of k increases at given a, the solutions 
will be progressively curved down both inward and outward 
of ro, to constitute a bundle tangent to the related powerlaw 
with slope given by Eq. (4). But, as stressed by Taylor & 
Navarro (2001), only a definite maximal value K cr i t can be 
withstood by the halo at equilibrium, lest the solution devel- 
ops pathologies like a central hole; this value marks a 'criti- 
cal' solution. Correspondingly, the values for K cr i t (a) may be 
found by a painstaking trial-and-error procedure, as done by 
Taylor & Navarro (2001) for the particular case a = 5/4. More 
conveniently, they may be computed directly through Eqs. (6) 



and (7); the former also provides the following exhaustive de- 
scription of the solutions. 

In the inner region only three situations possibly arise 
(Williams et al. 2004, Dehnen & McLaughlin 2005): (i) 7 
attains the value 70 asymptotically for r <C ro, which is un- 
physical since such a solution would end there with a diverg- 
ing gravitational force; (ii) the solution attains the value j a at 
a finite radius, and then continues to flatten down to the point 
of developing a central hole, which is obviously to be avoided; 
(iii) 7 approaches the value j a asymptotically for r <C ro, that 
is the only viable possibility. Thus physical solutions must 
approach the inner slope j a , and must steepen monotonically 
outwards. 

However, not all of them are tenable in view of their outer 
behavior (Austin et al. 2005, Dehnen & McLaughlin 2005). 
Here, again three situations formally arise that depend on the 
value of a: (i) 7 approaches the value 70 asymptotically for 
r ^> ro, which occurs for a > 35/27 and is unphysical since 
with such a solution the mass would diverge (in practice, it 
would be ill-defined, see § 5); (ii) 7 attains the slope 7/, at a 
finite radius and then goes to an outer cutoff, that occurs for 
a < 35/27 and is physically viable since the cutoffs actually 
set in at large radii exterior to the halo's bulk, as discussed 
below; (iii) 7 approaches the slope ji, asymptotically for r — > 
00 as it occurs for a = 35/27, which also provides a physically 
viable solution. 

2.2. Viable full solutions 

Summing up, the critical solutions satisfy both sets of cri- 
teria, but exist only for a < 35 /27 = 1 .296 (which excludes 
the case corresponding to the isothermal powerlaw); here- 
after these will be named ' a-profiles\ We compute them by 
integrating the second of Eqs. (6) with boundary conditions 
7 — > 7 fl and 7' — > for r — > 0. The results in the plane 7-7' 
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FIG. 4. — Characteristic radii (in units of ro) for the a-profiles as a function of a: radius r p where the circular velocity peaks, radius r,„ where the velocity 
dispersion peaks, radius where the slope of the density profile equals 3, and radius r;, where the slope of the density profile reads -74, =3(1 +a)/2. The blue 
dots refer to the solution with a = 5/4, and the red dots to a = 35/27. In the bottom right panel the red arrow indicates that the solution with a = 35/27 attains the 
slope 7j only asymptotically in a powerlaw tail. Note that as a increases, the peak of the velocity dispersion displaces relatively more than the circular velocity's. 



for various a between 1.2 and 35/27 are plotted as different 
lines in Fig. 1 (left bottom panel); there we can read out the 
values of 70 at r = ro (squares in the figure), and find the re- 
lated K cr j t after Eq. (7). The relation K cr it(aO is plotted in Fig. 1 
(right bottom panel); the resulting correlation is fitted by the 
expression K cr it(aO — 8.228-4.443 a to better than 1%. 

Our solutions include the prototypical one corresponding to 
a = 5/4 (highlighted in blue), first found by Taylor & Navarro 
(2001) and shown to be well fitted, if only in a region around 
ro, by the NFW formula (Navarro et al. 1997); for such a 
solution the precise value of K cr i t reads 2.674, see Dehnen 
& McLaughlin (2005). Our solutions also include the one 
corresponding to the limiting value a = 35/27 (highlighted 
in red) studied by Dehnen & McLaughlin (2005); for this 
K cl -it = 200/81 is found. 

In Fig. 2 we illustrate the radial dependencies of vari- 
ous interesting quantities for the a-profiles: density p and 
its logarithmic slope 7, circular velocity v\, velocity disper- 
sion a 2 , entropy K, mass M(< r), and gravitational potential 
$ = -J r co dr'GM(< r')/r' 2 ; the profil es are normalized to 1 
at the point ro within the halo's middle. Outwards of this, so- 
lutions corresponding to larger values of a have larger overall 
masses, flatter density profiles, higher velocity dispersions, 
and deeper potential wells. 

As mentioned above, the a-profiles with a < 35 /27 show 
an outer cutoff, to which the following remarks apply. First, 
the cutoffs generally occur beyond the virial radius (see § 5). 
Second, the velocity dispersion a 2 oc p 2 ^ K oc r a ~ 2j / 3 f a jj s 
rapidly outside of ro and vanishes at the density cutoff, thus 
preventing any finite outflow related to residual random mo- 
tions of DM particles. Third, the microscopic distribution 
function as obtained through the Eddington formula (e.g., 
Binney & Tremaine 1987) from pairs of p(r) and $(r) is re- 



alistic being nowhere negative (see also Zhao 1996; Widrow 

2000) . 

In Fig. 3 we show the ratio v 2 (r)/a 2 (r) entering Eq. (2); 
we also show the energy distributions associated to the 
a-profiles: the specific potential energy E pot = 4tt L dr' 
r' 2 p(r')v 2 (r')/M(< r), the random kinetic energy E^ n = 
4tt § Q r <\r'r l2 p(r') 3a 2 (r')/2M(< r), and the local virial ratio 
2£ , j c in/|£pot|- Note that the latter tends to 1 at large r, as is ex- 
pected for stable self-gravitating halos (see Lokas & Mamon 

2001) ; in fact, it closely approaches unity already for r w ro. 
In Fig. 4 we plot as a function of a various characteristic 

radii (in units of ro) of the a-profiles; specifically, we show 
the radius r p where the circular velocity v c (r) peaks, the radius 
r m where the velocity dispersion <r 2 (r) peaks, the radius rj 
where the slope of the density profile reads 7 = 3, and the 
radius r/, where the slope of the profile takes on the value 7& = 
3 ( 1 + a)/2. Note that r/, is always large, and formally recedes 
to infinity (as represented in the plot with an arrow) for the 
limiting case a = 35/27 that has a powerlaw falloff. 

For the sake of clarity, in these Figures we have confined 
the values of a to between 1.2 and 35/27; although the upper 
bound is already required for having no central hole and a fi- 
nite mass as discussed above, limited information existed so 
far concerning lower bounds from the isotropic Jeans equa- 
tion. 

3. THE ENTROPY SLOPES FROM 
COSMOGONIC BUILDUP 

Here we show that the cosmogonic histories of DM halos 
enforce additional, stringent limitations to the range of a; on 
the high side the resulting bound turns out to be pleasingly 
close to the value 35/27 derived above from the analysis of 
the static Jeans equation, but on the low side it constitutes a 
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new stringent result. 

Consider a DM halo with current bounding radius R, den- 
sity p, mass M oc pR 3 , velocity scales a 2 oc v 2 oc M/R, and 
entropy K cx a 1 / ' p 2 ! 3 . During the hierarchical buildup the run- 
ning growth rate scales as M cx pvR 2 with the inflow speed v 
set by the depth of the DM well and hence proportional to v 2 ; 
thus the above quantities may be expressed in terms of M and 
M to provide the scaling laws: 



RocM/M 2/3 



v 2 (xa 2 (y:M 2/3 



pozM 2 /M 2 



K(xM 4/3 /M 2/3 
(8) 



3.1. A heuristic approach 

Primarily we look for a link of the entropy slope a with the 
mass growth rate in simple terms; this analysis will have only 
heuristic scope, introductory to the full study performed in 
§ 3.2. Toward our scope, we make contact with the classic line 
of developments (Fillmore & Goldreich 1984; Lu et al. 2006) 
that analytically relate the collapse histories to the shape and 
growth rate of the primordial DM perturbations. The former 
may be expressed in terms of 5M/M cx M" e ; but since a halo 
virializes when 8M/M attains a critical threshold (e.g., 1.686 
in the critical universe), the shape parameter e = -dlog(l + 
z)/dlogM effectively governs also the (inverse) mass growth 
rate defined in terms of the scaling M oc (1 + z)~ 1 < /e . 

Then we use the standard relations between redshift z and 
cosmic time f, that may be expressed in the form (1 + z) oc 
t~ q . In the critical cosmology q = 2/3 holds; in the canonical 
Concordance Cosmology the value q m 2/3 still applies to 
z > 1, but as z decreases q progressively rises to values q w 0.8 
for z < 0.5. Then the mass growth history writes simply M oc 
t q l f - and the growth rate reads M/M = q/et; from Eqs. (8), 
simple algebra yields the scaling laws R oc t ( 2+ ?/ c V 3 ? p oc t~ 2 , 
and K oc f 2 < 1+ ?/ c )/ 3 . On eliminating t one obtains K oc R a and 
p oc /T 7 with slopes 



1 + 



1 



l+2e/q 



and 7 = 



6e 
2e + q 



(9) 



the latter being related to the former by the simple scaling 
link 7 = 6- 3a, as expected at r s=s ro from Eq. (4). From 
Eq. (9) a is seen to be slightly larger than unity by a quantity 
slowly depending one/q; this is consistent with K oc RM l l 3 oc 
R 2 ~~*/ 3 , the straightforward message from Eqs. (8) with M oc 

^3-7 

These relationships for a and 7 embody in a simple form 
some relevant cosmogonic and cosmological information. For 
Z > 1 where q w 2/3 applies, one explicitly has a w (2 + 
3e)/(l+3e) and relatedly 7 w 9e/(l + 3e), the latter being 
the well-known expression for the similarity solutions to the 
spherical infall model in a critical Universe. For z < 0.5 in- 
stead when q w 4/5 applies, one has a = (4 + 5e)/(2 + 5e) 
and 7 w 15e/(2 + 5e). As said, the remaining parameter e 
governs both the shape of the initial perturbations and the 
times of their collapse. For instance, if the standard top-hat 
shape applied with e s=s 1 corresponding to a monolithic col- 
lapse, from the above relations one would obtain a « 1 .25 and 
7 w 2.25 for halos that virialize at z > 1, as pointed out early 
on by Bertschinger (1985); for more massive halos that virial- 
ize at z < 0.5, from the scaling laws one anticipates a slightly 
steeper a w 1.29 and a. flatter 7 w 2.15. 

Qualitatively, different values of the slope 7 may be related 
to other values of e, such as occurring in an articulated shape 



of the initial perturbations, that also yields a sequence of col- 
lapse times. With a realistic bell-shaped perturbation e s=s 1 
still applies to its middle; in the outer regions e > 1 implies 
the accretion to become slower, and correspondingly the den- 
sity slope 7 to steepen there relative to the middle one. The 
central regions are described by e <C 1, corresponding to fast 
collapse; the scaling law then suggests the density slope to be 
flatter there. 

However, this simple approach cannot be carried very far. 
At small radii also angular momentum effects come into play 
to lower the density profile; they introduce another radial scal- 
ing related to centrifugal barriers, as discussed in detail by 
Lu et al. (2006); these authors show how a density slope 
around -1 is enforced if efficient isotropization of particles' 
orbits takes place during the fast collapse, while still flat- 
ter slopes occur with more even angular momentum distri- 
butions. These, however, must be consistent with the overall 
equilibrium requirements; these orbital features are actually 
subsumed into the macroscopic Jeans Eq. (2) with K oc r a 
and provide the running 7(7-) = 3 a/5 + 3 v 2 (r)/5 a 2 (r) (the last 
term is shown in Fig. 3) that is to replace the stiff scaling value 
7 = 6-3 a above, but still recovers it at the single radius ro (cf. 
Eqs. [4] and [6]). 

So we will depart from the above heuristic approach, and 
focus on computing the values of a to be inserted into Jeans, 
also implementing precise mass and time dependencies of the 
growth rate e. To focus these sensitive issues a refined analysis 
and a careful discussion are required and performed next. 

3.2. Advanced, semianalytic analysis 
We derive a robust time evolution and mass dependence of 
M from an updated version of the Extended Press & Schechter 
formalism (EPS, see Lacey & Cole 1993). 

We are mainly interested in the average evolutionary be- 
havior of halos, so we compute the mean growth rate M of the 
current mass M as 



/•M 

M = dM' (M-M') 
Jo 



A 2 P M 



AM' At 



(10) 



in terms of the differential merger rate & 2 Pm>^,m /AM' At 
given in Appendix A. This incorporates the full cosmological 
framework, and the detailed cold DM perturbation spectrum; 
as to the former we implement the Concordance Cosmology, 
and as for the latter we adopt the shape given by Bardeen et al. 
(1986) corrected for baryons (Sugiyama 1995), and normal- 
ized as to yield a mass variance <7& = 0.8 on a scale of 8/z -1 
Mpc. We improve over previous studies of this kind (e.g., 
Miller et al. 2006; Li et al. 2007; Salvador-Sole et al. 2007) 
by using the merger rate appropriate for ellipsoidal collapses 
(Sheth & Tormen 2002) as provided by Zhang et al. (2008); 
this is shown by the latter authors to yield formation times 
and growth rates in much closer agreement with numerical 
simulations than the standard EPS theory. Further details are 
presented in Appendix A. 

Eq. (10) provides the relation M(M,t) illustrated in Fig. 5. 
This may be viewed as a first order differential equation for 
M(t), that is numerically solved once the mass M(to) at the 
present time to has been assigned in the way of a boundary 
condition. In terms of M(t) so computed, the evolution of 
the halo properties are obtained through Eqs. (8); the results 
are illustrated in Fig. 6, where we plot the time and redshift 
evolutions for the inverse growth rate e, the overall radius R, 
the mass M, the characteristic velocities v oc a, the entropy 



DM EQUILIBRIA IN GALAXIES AND GALAXY SYSTEMS 



7 




FIG. 5. — The relation M(M,z) from the average cosmogonic evolution, see § 3.2 for details. 



K, and finally the corresponding slope a = dlogTT/dlog/?. 
The curves are for different current values of the halo masses 
M(t ) ranging from 10 11 to 10 15 M Q . 

In fact, in these detailed mass accretion histories one can 
recognize a transition epoch when the circular velocity v 2 of 
the halo attains its highest value along an accretion history 
(see Li et al. 2007). The transition separates two stages: an 
early stage of fast collapse, when e < 0.4 corresponds to a 
growth rate still high on the running Hubble timescale, and 
nearly constant; and a later phase of slower accretion when e 
takes off from the value 0.4, i.e., the growth rate is low and 
decreasing. Small halos with M < 10 12 M Q have high tran- 
sition redshifts, hence are now well advanced into their slow 
accretion stage; on the other hand, very massive ones with 
M > 1O 15 M have low transition redshifts, so today they have 
just completed their fast accretion stage. A similar behavior 
has been identified in a number of recent numerical simula- 
tions (Zhao et al. 2003; Diemand et al. 2007), and has been 
recognized also in other semianalytic computations based on 
the simpler, less precise EPS theory (see Li et al. 2007; Neis- 
tein et al. 2006). 

3.3. Results 

From Fig. 6 we also see that the values of a at the tran- 
sition are narrowly constrained to within 1.27 ±0.02, i.e., a 
1.5% range. We note that this range explains why narrowly 
constrained phenomenological values of a are found within 
the halo bulk from many and different simulations since the 
year 2000. In particular, our values turn out to agree with 
the best fitting exponent of Q(r) [or K(r)] with its scatter re- 
cently given in detail by Ascasibar & Gottlober (2008) within 
the halo bulk during calm infall; these authors warn against 
extrapolating their results to early times and outer radii. 

Correspondingly, in our findings the values of a before the 
transition should not be overinterpreted nor inserted into the 
Jeans equation. In fact, at these early epochs the halos expe- 
rience rapid changes and are clearly out of equilibrium; dur- 



ing this stage the N-body experiments (see Zhao et al. 2003; 
Peirani et al. 2006; Hoffman et al. 2007) highlight turmoil 
associated with fast collapse and major mergers, effective dy- 
namical relaxation and mixing of DM particles. Taking up 
from § 3.1, it is during these events that the orbits' angular 
momentum must set to a near isotropic distribution such as to 
lower the inner densities. 

But as the turmoil associated with the fast collapse subsides, 
the orbital structure is bypassed and replaced by the macro- 
scopic entropy distribution K oc r a with uniform values of a, 
as phenomenologically found to hold throughout the now sta- 
ble halos' bulk. Then the values of a at the transition epoch 
and the corresponding edge R = ro can be derived from the 
scaling laws and applied to the whole inner region, to obtain 
7(r) = 3 a/5 + 3 v 2 (r) / 5 a 2 (r), taking on the value 6 - 3 a only 
at the point ro according to the results of § 2. 1 . 

The value of a is expected to persist in the bulk during the 
subsequent development in the slow accretion stage, when 
mass is added smoothly and mainly onto the outskirts, with 
little mixing. Such a development is borne out by the widely 
shared scenario of halo buildup from the inside out, champi- 
oned by Bertschinger (1985), Taylor & Navarro (2001), and 
many others up to Salvador-Sole et al. (2007). Specifically, 
the analysis by the former authors concerning the randomiza- 
tion and the stratification of particles' orbits during secondary 
infall indicates that their internal disposition is little affected 
by subsequent, smoother mass additions on top of a closely 
static potential well. Thus during the late, slow accretion stage 
the structure for r < r will be unaffected. 

Meanwhile, the outer scale R oc M/M 2 / 3 beyond ro is in- 
creasingly stretched out to substantially larger values (see 
Fig. 6, top right panel); correspondingly, we find the outer 
slope a(R) = dlogTT/dlog/? at R > ro to decline in time 
(see Fig. 6, bottom right panel). This translates into a ra- 
dial decrease a(r) on maintaining entropy stratification also 
into these regions (see Bertschinger 1985, Taylor & Navarro 
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FIG. 6. — Evolution of the various quantities along the average accretion history of DM halos: mass M, radius R, velocity v, (inverse) mass growth rate 
e = — dlog(l +z)/dlogM, entropy K, and slope a = dlog K/dlogR. The curves are for different current halo masses in the range 10" - 10 15 M©; from top to 
bottom, blue lines refers to 10" M©, cyan to 10 12 M©, green to 10 13 M©, orange to 10 14 M©, and red to 10 15 Mq. Lines are dashed in the fast accretion regime, 
and dotted in the slow accretion phase; the solid line marks values at the transition epoch, see § 3.3 for details. In the bottom right panel, the shaded area 
highlights the narrow range expected for the values of a at the transition. 



2001). Such a decrease is to be expected also from consider- 
ing that K(r) oc r 2 " 7- / 3 — > r applies, on the basis of the related 
asymptotic behaviors of p(r) oc r~ 7 and er 2 (r) oc r 2 " 7 where 
7^3. We find the values of a to decrease by less than 2% 
out to r w 2 ro, and less than 10% out to r w 5 ro. 

Concerning the Jeans equation, this too (with its boundary 
conditions set at r = 0, see § 2.2) technically works from the 
inside out, and may be extended to the range r > ro with de- 
creasing slope a(r). The corresponding density run is some- 
what steepened for r > ro, as if shifting from the upper to the 
lower curves in Fig. 2; we have checked that the outer density 
profile lowers by less than 20% for r < 5 ro. State-of-the-art 
simulations still yield little information on the outskirts for 
r > ro owing to low densities and noisy data there discussed 
by Ascasibar & Gottlober (2008). We present elsewhere the 
complete solutions outside ro as a testbed for future simula- 
tions aimed at resolving the outskirts, and at fitting gravita- 
tional lensing observations (see discussion in § 5). 

In a similar vein, we predict a minor increase of a at the 



transition from 1 .25 to 1 .29 in moving from galaxy- to cluster- 
sized current masses. A fit accurate to better than 0. 1 % to this 
average relation is provided by a ~ 2.71-0.23 log(M/M©) + 
0.009 [log(M/M Q )] 2 ; in particular for galaxy clusters, this im- 
plies values of a » 1 .27 for Abell richness classes r w and 
a w 1.29 for richness class r > 3. The correlation a(M) may 
help to understand the differences in the values of a found 
in numerical simulations; in fact, the highest value 1.29 to 
date has been reported by Rasia et al. (2004) in looking at 
cluster-sized halos, while the value 1 .25 first found by Taylor 
& Navarro (2001) refers to galaxy-sized halos. Observational 
evidence supporting the above picture is presented in § 5. 

To sum up, from the cosmological buildup we find for the 
entropy slope the narrow allowed range 

aw 1.27 ±0.02; (11) 
this we use in our Jeans-like description of DM equilibria 
within the halos' bulk. The bounds on a are relevant in two 
respects: the upper one strengthens the formal limit 35/27 = 
1 .296 required for the solutions of the static Jeans equation to 
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have viable profiles in the far outskirts (see § 2.2 and discus- 
sion in § 5); the lower one complementarily restricts the range 
of the viable values. 

4. A PHYSICAL MEANING FOR K crit 

Here we propose to discuss the physical meaning of the 
quantity K cr [ t (defined in § 2.1 following Taylor & Navarro 
2001) in terms of a direct energetic argument. We first note 
that given the asymptotics of the a-profiles, the velocity dis- 
persion v 2 (r) must peak at a point r p > ro (see Fig. 4) to the 
value v 2 = 4TrGp(r p )r 2 ; at this easily identified point the con- 
stant k also writes 4irGp(r p )r 2 /a 2 = v 2 p /a 2 , and thus is seen 
to compare the estimate a 2 for the random kinetic energy with 

the estimate v 2 for the gravitational potential. We then recall 
from § 2.1 that a halo in the viable equilibrium condition for 
a given a is marked by a maximal value K cr it(a), that is, by a 
minimal level of randomized energy a 2 = v 2 /K cr \ t . 

We now show that to a very good approximation the value 
for Kent may be obtained from considering incomplete con- 
version of the infall kinetic energy v? nf / 2 of DM particles into 
their dynamically randomized energy 3 a 2 / '2 during forma- 
tion; a process introduced by Bertschinger (1985), different 
from the complete violent relaxation proposed by Lynden- 
Bell (1967) and discussed widely, including the recent Arad 
& Lynden-Bell (2005), Trenti & Bertin (2005) and Trenti et 
al. (2005). 

We propose to derive the critical value of k = v 2 /a 2 from 
the relation 



3 <7„ 



with 



77 = 



v 2 P /vl 



inf 



2 l + v 2 p /v{ 



(12) 



inf 



This expression for the effective conversion efficiency r] has 
been constructed on considering that: (i) it has to depend on 
the ratio of the (1-D) infall velocity squared v 2 nf to the depth 
of the potential well as measured by the (2-D) circular veloc- 
ity squared v 2 ; (ii) the formal limit v 2 =2a 2 corresponding to 
k = 2 (the isothermal powerlaw, see §2.1) must be approached 
for a fast inflow with v? f ^ v 2 , which quantifies the incom- 
pleteness of the energy conversion during the transit of a bullet 
with crossing time short compared with the halo's dynamical 
time. 

As a simple example guided by the canonical theory of 
gravitational interactions (see Saslaw 1985; Cavaliere et al. 
1992) we may consider that vf nf /v 2 w 3 is to hold for having 
considerable energy transfer; using this simple approximation 
in Eq. (12) yields K a i t i=s 2.6, a value pleasingly close to 2.64 
as obtained numerically in § 2.2 for the particular a-profile 
with a = 1 .25. On the other hand, this simple argument cannot 
pinpoint values of K cr ; t to the precision required to discrimi- 
nate different values of a. 

So, we may test our proposed rj as given in Eq. (12) for 
different values of a by using instead energy conservation to 
express the infall velocity v 2 nf « 2 A(j> p , in terms of the adi- 



mensional potential drop A<j> pr (normalized to v 2 ) from the 
turnaround radius to r p (see Bertschinger 1985; Lapi et al. 
2005). Inserting this into Eq. (12) yields the expression 



v 2 

_ p 



-4 = 2+- 



1 



A(j) p 



(13) 



this highlights the deviation of k ci h from the value 2 that 
would apply for the isothermal powerlaw (see § 3). Using the 
runs of (f> appropriate for the a-profile as given in Fig. 2, we 
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FIG. 7. — The /t cr i[(o) relation, as in Fig. 1 (bottom right panel). Crosses 
illustrate the goodness of our approximation to re cr jt(») based on a direct en- 
ergetic argument, see § 4 for details. 



compute values for K cr ; t as a function of a; note that these de- 
crease for increasing a which implies larger A<fi pr , see Fig. 2. 
Fig. 7 shows how well Eq. (13) recovers the K CI \ t (a) relation 
obtained through numerical integration of the Jeans equation. 

This excellent fit to n all means that the expression for -q 
in Eq. (12), more than a successful guess, is describing what 
actually occurs in the macroscopic equilibria expressed by the 
Jeans equation in keeping with the numerical simulations. At 
the referee's suggestion we add that what ultimately the first 
Eq. (12) describes with the isotropic form of its l.h.s. is a 
partial effective conversion at r 
tangential motions. 



r p from radial infall into 



5. DISCUSSION AND CONCLUSIONS 

The general picture of DM halo formation as it emerges 
from recent Af-body simulations envisages a progressive 
buildup through a two-stage development. At early epochs 
major mergers associated with fast reshaping of the bulk grav- 
itational potential cause particles' orbits to dynamically re- 
lax and mix into a definite equilibrium disposition. At later 
epochs the halo outskirts progressively build up through rare 
if persisting minor mergers and by smooth mass accretion, 
settling to a quasi-static equilibrium. 

The stable equilibria soon after the early stage are effec- 
tively described in terms of a phenomenological entropy run 
K (xr a with uniform slope a somewhat larger than 1 . Spurred 
by this powerlaw behavior, one might look for powerlaw so- 
lutions of the Jeans equation for the equilibrium density pro- 
files, only to find these just tangent to the actual solutions. 
From here our analysis goes beyond powerlaw densities, and 
provides two related blocks of novel results. 

First, we have computed and illustrated a wide set of novel 
solutions to the isotropic Jeans equation, that we name a- 
profiles, including as two particular instances the prototypes 
investigated by Taylor & Navarro (2001) and Dehnen & 
McLaughlin (2005). We have stressed that at given a the so- 
lutions form a tangent bundle with p(r) progressively curved 
down as the gravitational pull measured by k = v p /cr 2 in- 
creases above the value corresponding to the pure powerlaw 
profile 7 = 6 - 3 a; from such bundles the a-profiles stand out 
for having a shape maximally curved to the point of comply- 
ing with sensible outer and inner behaviors, i.e., finite total 
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FIG. 8. — Analytic fits to the density runs (top panels) and residuals (bottom panels) for the a-profiles with a = 1.25 (left panels) and a = 35/27 (right panels). 
Solid lines illustrates the exact solutions from numerical integration of the Jeans equation, dashed lines show our analytic fitting formula Eq. (Bl), and dotted 
line (only in left panels) is the fitting formula by Taylor & Navarro (2001, their Eq. [10]); note that Eq. (Bl) for a = 35/27 coincides with the exact solution given 
by Dehnen & McLaughlin (2005). 



mass and no central hole. This occurs for a maximal value 
Kcrit(oO °f the gravitational pull v 2 relative to the random com- 
ponent <Tp, i.e., for a minimal a 2 compared to v 2 ; this also cor- 
responds to the maximal slope j p = 3 a/5 + K cr [ t at the point r p . 
Such solutions still constitute a large set, being parameterized 
in terms of a wide overall range 1 < a < 1 .296 for the values 
of a. 

Second, we have shown that toward further constraining 
these values the key step is provided by the halos' cosmolog- 
ical buildup; in fact, the scaling laws of § 3 narrowly restrict 
the viable range to a w 1.25 — 1.29 w 1.27 ±0.02, correspond- 
ing to current masses M <~ 10 11 - 10 15 M Q . We have found 
these bounds to apply at the epochs of transition from fast col- 
lapse to slow accretion, when the potential wells settle to their 
maximal depths along an accretion history (see § 3). In fact, 
our semianalytic analysis not only recovers the two-stage de- 
velopment observed in many recent A^-body simulations, but 
also explains why only such narrowly constrained values of 
a are actually found in many different numerical works since 
the year 2000. The values K cr ; t is correspondingly constrained 
to the range 2.7 > K cl i t > 2.5; we trace back these values to 
limited randomization for the infall kinetic energy of DM par- 
ticles during accretion, which leads to minimal values of a 2 . 

The a-profiles from the Jeans equilibrium with a and K cr i t 
so constrained imply that galaxies have a generally more com- 
pact central structure than groups and clusters, with apprecia- 
ble deviation from self-similarity; on the other hand, the for- 
mer have currently developed extensive outskirts, while the 
very rich clusters do not yet. 

To embody in a handy expression the basic features of the 
density runs for the a-profiles, we provide here an analytic 
fitting formula of the form often advocated on phenomenolog- 
ical grounds (see Zhao 1996; Kravtsov et al. 1998; Widrow 



2000; Barnes et al. 2006), that reads 

p(r) K TTl^TTTTTW' (14) 

(r/r Q ) [l + w (r/ro) \ 

with 7 fl governing the inner profile, q the outer one, w and u 
the position and sharpness of the radial transition. Here for 
the first time such a profile is substantiated with parameters 
expressed in terms of the derivatives of the Jeans equation 
(see § 2.2) and of the values for a as detailed in Appendix 
B. We stress that the Jeans equlibrium (see Fig. 2) requires 
that the central slopes are always flatter, and the outer ones 
always steeper than in the empirical NFW fit. The goodness 
of the above fitting formula is discussed in Appendix B, and 
illustrated in Fig. 8. 

We probe the physical relevance of the a-profiles by com- 
paring in Fig. 9 the related projected density runs with the 
results from recent, exquisite observations of strong and 
weak gravitational lensing in and around massive galaxy clus- 
ters, which just require flat inner, and steep outer behaviors 
(Broadhurst et al. 2008, Umetsu & Broadhurst 2008, and ref- 
erences therein). It is seen how straightforward are our fits in 
terms of high-mass a-profiles that naturally exhibit just such 
features, compared to those in terms of the NFW empirical 
formula that not only has problems of ill-defined mass, but 
also would require particularly high concentrations and old 
ages. Our findings substantiate the a-profiles, and stress their 
value as to specific predictions for future simulations aimed 
at resolving, and for weak gravitational lensing observations 
aimed at probing the outer cluster densities. A full discussion 
of these issues deserves a dedicated paper. 

Two comments are in order concerning the validity of the 
Jeans equation in the simple form of Eq. (2). First, reason- 
able anisotropics in the DM are described by the standard 
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FIG. 9. — Surface density runs from the ^-profiles for a = 1.25 (blue), 
a = 1.27 (green) and a = 1.29 (red) and from the empirical NFW profile 
with concentration parameter c = 5 (black dotted), compared with a fit to the 
gravitational lensing data by Umetsu & Broadhurst (2008, black dashed) in 
terms of an NFW profile forced to high concentrations c as 14, but actually 
overshooting the central data. Note how our physical profiles do better both in 
the central and in the outer region while relieving, along with the two-stage 
development, the problems of high concentrations and short ages that arise 
from NFW fits as discussed by Broadhurst et al. (2008). 

parameter = 1 — Og/of (Binney 1978), which numerical 
simulations suggest to increase with r from central values 
0(0) « -0.1 -j- meaning near isotropy to outer values close 
to 0.5, meaning progressive outward prevalence of radial mo- 
tions (see Hansen & Moore 2006). Such anisotropies inserted 
into Jeans produce small variations of the a-profiles, to the 
effect of slightly increasing the limiting value a = 35/21 - 
4 /3(0) /27; thus in the related limiting solutions the inner slope 
is flattened 7fl = 3 a/5 + 6 0(0) /5 = 7/9 + 10/3(0) /9, the mid- 
dle slope 7o = 6-3 a = 19/9+4/3(0)/9 is steepened, and the 
outer slope 7/, = 3 ( 1 + a) /2 = 3 1 /9 - 2 /?(0)/9 is also steepened 
with respect to the isotropic case, see Dehnen & McLaughlin 
(2005). The oppositely extreme case of total radial anisotropy 
with constant = 0(0) 1 discussed by Austin et al. (2005) 
would yield a lower limiting a = 31/27 ~ 1.148, while also 
producing a steep inner slope j a = 17/9 w 1.88 (cf. § 2.1). In 
brief, the distinctive features of the a-profiles are somewhat 
enhanced, if anything, with radially increasing anisotropies. 

Second, the static Jeans equation becomes increasingly vul- 
nerable in going into the outermost regions, since static or 
even dynamic equilibrium is not granted beyond the virial ra- 
dius. Specifically, in the Concordance Cosmology the Hubble 
expansion takes over at radii exceeding a few times the virial 
radius where 17a (#o>" /v c ) 2 > 1 holds. This weakens the case 
for the formal upper limit a = 35/27 as derived from Jeans 
analysis but correspondingly highlights the relevance of the 
ceiling a w 1.3 for a few 1O 15 M halo we have derived from 
cosmogonic evolution. On the other hand, any halo termi- 
nation will be uncertain by full factors of 2; with flat outer 



density slopes, as the case with a > 35/27 (see end of § 2.1), 
the related total masses will be ill-defined by similar factors, 
corresponding for galaxy clusters to a disturbing uncertainty 
by ±2 in Abell richness class (see Yee & Lopez-Cruz 1999). 

Finally, we comment on our use of the DM 'entropy' K(r) = 
cr 2 (r)/ p 2 ^(r) - actually the analog of the thermodynamic adi- 
abat, related to entropy by s = 3/2 log K - rather than the 
equivalent functional 'phase-space density' Q = K~ 3 I 2 - that 
underlies the 'information' log Q = —s with the obvious mi- 
nus sign. Matter-of-factly, the pressure gradient d(pa 2 )/dr 
that withstands the self-gravity in Eq. (1) is most directly ex- 
pressed in terms of the macroscopic quantity K, in the form 
d (p 5 / 3 K)/dr. In the way of interpretations, our lead is pro- 
vided by the positive-definite gradient As = 3/2 AlogA', and 
specifically by its adimensional version dlog/sT/dlogr = a. 

We find As to be convenient in understanding the empiri- 
cal evidence as to positive and uniform of a within the halos' 
bulk after the transition; then and there, in fact, the particle 
orbits settle to a definite dynamical regime, relaxed toward a 
periodic, ordered behavior (Bertschinger 1985), correspond- 
ing to lower entropy layers at the center. In the outer regions, 
on the other hand, the particle orbits stratify but mix poorly 
during the late, slow buildup (see Bertschinger 1985; Tay- 
lor & Navarro 2001) with their apocenters distributed over 
larger and larger sizes; we understand this dynamical regime 
in terms of K or As spread over increasingly stretched radii 
R > r p , resulting in flatter slopes. 

It rests with the numerical values of a to translate the 
fine orbital details into the macroscopic average of Jeans, as 
clearly indicated by the consistency of the results from the 
latter with semianalytic cosmogony and with simulations. 

We end by stressing that the formulation of DM distribu- 
tion in terms of 'entropy' - despite the standard objections 
about defining entropy for a collisionless medium out of ther- 
mal equilibrium and dominated by long range weak gravity 
- opens an interesting perspective; this envisages direct com- 
parisons (see Faltenbacher et al. 2007) with the truly ther- 
modynamic entropy of the other major component in galaxy 
systems, namely, the collisional plasma in local thermal equi- 
librium constituting the intracluster medium. In fact, in the 
outer regions these entropies though related as being both 
originated by conversion of the infall energy under weak grav- 
ity, differ owing to their origin from different conversion pro- 
cesses (Lapi et al. 2005; Cavaliere & Lapi 2006); in addition, 
in the central region of most clusters the hot plasma is heated 
up by the activity of supermassive black holes originating un- 
der strong gravity conditions. 
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APPENDIX 

A. MERGING RATE OF DM HALO UNDER ELLIPSOIDAL COLLAPSE 

To reliably compute the growth rate M in Eq. (10), we use the merging rate of DM halos under ellipsoidal collapse as provided 
by Zhang et al. (2008); this reads 



d 2 P, 



aoe 



-a\Aa-/2<j- M 



dM'dt 



/2^(A(j2)3/2 



da 2 



dM' 



d5 r 



dt 



1 +ai 



Aa 2 



' M 



3/2 



1+: 



«1 



'Act 2 



T(3/2) 



M 



(Al) 



12 LAPI & CAVALIERE 

where Act 2 = a 2 M i - ct^ is the difference of the mass variances corresponding to the masses M' and M, and 

ao = 0.866 {l-0.133 [<5 c (z) 2 /^m]"°' 615 } , 

a x = 0.308 [5 c (z) 2 /a 2 M ]^ U5 , (A2) 
fl2 = 0.0373 [^(zf/CT 2 ,]" ' 1 ' 5 . 

Note that the standard EPS expression for the spherical collapse model (Lacey & Cole 1993; Kitayama & Suto 1996) is recovered 
on setting ao = 1 and a\ = a-i = 0. 
Eqs. (Al) and (A2) involve the critical threshold S c (z) for collapse, which is computed according to 

*c(z) = <5 c ^ (A3) 
D(z) 

in terms of the normalization value S c w 1.686 [1 +0.0123 logO^z)] and of the growth factor 



2 1+z 



4 + T^^ (z) -T^o - (z)+0 - 7(z) 



(A4) 



in the above the quantity Qm(z) = (1 +z) 3 /[ 1 - ( 1 +z) 3 ] represents the evolved matter density parameter of the Universe 

in the Concordance Cosmology. 
Eqs. (Al) and (A2) also involve the mass variance of the primordial perturbation spectrum, that writes 

<j 2 m = Ait / dkk 2 P(k) , (A5) 
Jo 

where k s = (6tt 2 p M / 'M) 1 / 3 is the wavenumber corresponding to the mass M, and p M ~ 2.8 x tt M h 2 M Q Mpc~ 3 is the mean 
background matter density of the Universe. The mass variance is then normalized to the value <j% = 0.8 for the mass corresponding 
to a scale of 8 h~ l Mpc. As for the power spectrum, we adopt the shape 



P(k)(xkT 2 (k) (A6) 



in terms of the transfer function 



T(k) ^° g(1 2 + 3 2 4 3 q 4q) [l + 3.89< ? + (16.1 9 ) 2 + (5.46 9 ) 3 + (6.71 ? ) 4 ]- 1/4 (A7) 

given by Bardeen et al. (1986). Here q = k/{hT) includes the shape correction factor Y = Q M he~ nh ~^ 2 '* < - ni '/ nM ^ due to baryons, 
after Sugiyama (1995); we recall that the values h = 0.72, Q M = 0.27, and Qj = 0.044 are adopted throughout the paper. 

B. FITTING THE DENSITY RUNS OF THE a-PROFILES 

Here we provide an explicit analytic fit to the a-profiles, whose numerical computation is sensitive to precise initialization. We 
use the expression 

p(r)ozx-~<"(l + wx u y q (Bl) 

inspired by the exact solution for a = 35/27 given by Dehnen & McLaughlin (2005), and representing an extension of the NFW 
(Navarro et al. 1997) formula; in the above x = r/ro is the adimensional radius, the inner slope has been fixed to j a = 3 a/5, 
while m, q, and w are parameters dependent on a to be determined. 

To this purpose, we impose that at r = ro the density slope from Eq. (Bl) and its first and second logarithmic derivatives equal 
those of the a-profiles. Thus we find 



w u 2 qw , m 3 qw(l-w) 

la + U q—= lQ 77^2=70 * , 3 =7o • (B2) 

w+1 (1 + wY (l + wy 

Then we solve the above expressions for w, u, and q, to obtain 



, 7 (70-7a) (70 -la) (70-7a)(w+l) 

w = 1 - -2 n q = ; — u = . (B3) 

7o w % 1 W 

The values of j' and j{j are easily computed as a function of a. First, Eqs. (6) evaluated at r = r give 

3 2 

70 = 5 ^crit(a) + - (70 " la) (70 - lb) 



(B4) 



2 

7o = 370(270 ~ia- ib) ; 
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then these quantities are easily recast in terms of a on recalling from the main text (§ 2.2) that j a = 3 a/5, 70 = 6- 3 a, 72, = 
3(1 + q)/2, and K crit (a) ~ 8.228-4.443a. 

To sum up, Eqs. (B3) and (B4) provide the parameters w(a), w(a), and <7(a) as a function of a. E.g., for a = 35/27 we obtain 
w = 1 .00, u = 0.44, and q = 6.00; thus our fitting formula recovers the exact solution by Dehnen & McLaughlin (2005). At the 
other extreme, for a = 1.25 we obtain w = 0.23, u = 0.39, and q = 20.77; here our fitting formula approximates well (see Fig. 8) 
the exact solution, improving upon the fit by Taylor & Navarro (2001, their Eq. [10]) 4 . In fact, we find that for any a between 
1.25 and 35/27 our fit works to better than a few % throughout the significant radial range from 10 _1 to several ro; note that 
both the radial range and the quality of the fit increase for increasing a, up to the point of pinning down the exact solution for 
a = 35/27. 
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